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We continue the discussion of the fermion models on graphs that started in the first paper of 
the series. Here we introduce a Graphical Gauge Model (GGM) and show that : (a) it can be 
stated as an average/sum of a determinant defined on the graph over Z2 (binary) gauge field; (b) 
it is equivalent to the Monomer-Dimer (MD) model on the graph; (c) the partition function of the 
model allows an explicit expression in terms of a series over disjoint directed cycles, where each 
term is a product of local contributions along the cycle and the determinant of a matrix defined 
on the remainder of the graph (excluding the cycle) . We also establish a relation between the MD 
model on the graph and the determinant series, discussed in the first paper, however, considered 
using simple non-Belief-Propagation choice of the gauge. We conclude with a discussion of possible 
analytic and algorithmic consequences of these results, as well as related questions and challenges. 

PACS numbers: 02.50.Tt, 64.60.Cn, 05.50. -|-q 

Gauge theories, stated in terms of fermions and gauge fields (e.g. associated with a vector potential), are common in 
theoretical and mathematical physics jlj|2j- Normally in physics, e.g. discussing quantum electrodynamics or quantum 
gravity, these popular theories are defined over continuous spaces, or their natural discretizations, e.g. triangulated 
Eucledian grids. In the discretized versions, e.g. in Lattice Gauge Theories [S , fermions are normally associated with 
vertexes of the grid, while gauge variables reside on edges. 

In this paper we extend this standard discretized construction to arbitrary graphs and show that the gauge theory 
approach, native of physics, can be useful for getting new nontrivial relations between different graphical models that 
describe computer science problems defined on arbitrary graphs. We introduce and discuss a Graphical Gauge Model 
(GGM). Gauge fields in our construction correspond to standard binary variables, which could also be called Ising 
spins/variables, or more formally, the gauge group of the theory is Z2. Two objects emerging in any gauge theory, 
determinants and loops, are therefore natural participants of our description. We also find that this approach and 
language fits naturally with the Loop Calculus introduced in [T, '5] and extended in the first paper of the series fB] to 
the Gaussian Graphical models on graphs. 

The power of GGM is in its natural operational flexibility: changing the order of integrations and modifying the 
integrand in the expression for the model's partition function result in a variety of nontrivial relations, some of them 
discussed in this paper. Integration over the Grassman-fermion variables turns the partition function of GGM into a 
(■-function dependent on the gauge (binary variable) configuration. Here the ^-function is understood as a generating 
function for the expectation values of the Grassman variables and their combinations. In this formulation it is related 
to the inverse of the Ihara (^-function of the graph [7, . 

Even though we are making a point in promoting the language of Grassman/fermion integration in this paper, and 
the series in general, our two main statements are made in terms of "normal" objects, e.g. determinants, disjoint 
oriented cycles, and also partition functions of the monomer-dimer model. Given that this latter object did not 
appear in the first paper of the series, we find it useful to state it casually right away (see the rhs of Eq. (|6| for 
formal definition). Consider a graph, and cost functions, Wa and Wab, associated with the vertices and edges of the 
graph. A monomer-dimer configuration on the graph is a set of colorings of vertices and edges so that either any 
vertex of the graph is colored and then no adjusted edges are colored, or the vertex is not colored but then one of the 
adjusted edges is colored. The partition function of the monomer-dimer model on the graph is the sum over all allowed 
monomer-dimer configurations/colorings, where each individual contribution is a product of factors associated with 
colored vertices and edges over the graph ^. Armed with this definition let us now state the main results reported in 
the paper: 

• The partition function, Zmjj, of the monomer-dimer model on a graph Q is expressed in terms of a matrix 
H built from the monomer and dimer weights placed at the diagonal and off-diagonal elements, respectively. 
Specifically, Zmd{Q) is stated as a series over the oriented disjoint cycles C G OT)C{Q) of the graph. An oriented 
disjoint cycle C G OT>G{Q) is represented by a disjoint union of simple oriented loops. Each term in the series is 



Note, that "dimers" and "monomers" are terms used in statistical physics which are also fully equivalent to "perfect matchings" and 
"imperfect matchings" in the terminology commonly accepted in computer science, see e.g. [8]. 
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FIG. 1: Schematic description of relations between models and concepts discussed in the paper. 



equal to the determinant of the original matrix H with the cycles excluded, H\ 



multiplied by the product 



along the simple loops of the cycle of the corresponding off-diagonal elements taken with the reversed signs. 
(See Eqs. ( |6|l3| .) 



• The determinant of H is stated as a series over oriented disjoint cycles C E ODC{Q) of the graph Q, where each 
term is equal to the partition function, Zmd{Q \ of the monomer-dimer of the original graph with the cycle 
excluded, multiplied by the product along the simple loops of the cycle of the off-diagonal elements taken with 
the reversed signs. (See Eqs. pSje] ).) 



Three remarks are in order. First, the two main statements are equivalent, in fact, one is a kind of an inverse of 
the other (see Appendix |A] for clarifications). Second, the first statement has an immediate algorithmic consequence: 
it may be used for an approximate computation of the monomer-dimer partition function (which is known to be an 
T^P-complete, i.e. a counting problem of the likely exponential complexity [3]) via a truncation of the determinant 
series (the complexity of a determinant evaluation is cubic in the graph size). Third, the statement number two 
(expansion of the determinant in a series over oriented disjoint cycles) can be derived by implementing the Gauge 
fixing approach in the spirit of [H [5], however selecting a gauge different from the Belief Propagation (BP) gauge. 
The latter resulted in the Loop Series expansion for the determinant, described in the first paper of the series [B]. 

A schematic set of relations between the two main statements and other results and models discussed in the paper, 
are shown in Fig. [l] The distribution of material is as follows. GGM is introduced in Section TK Direct relations 
between GGM and partition function of MD and DCS over determinants are established in Section |IB| and Section 
|IC| respectively. Inverse of the relation, expressing determinant as a series over partition functions of MD models 
on the original graph with disjoint cycles excluded is discussed in Section III] with some auxiliary material placed in 



Appendix Kl and Appendix p| Section III is reserved for Summary and Conclusions. 
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I. GRAPHICAL GAUGE MODEL, MONOMER-DIMER AND DISJOINT CYCLE SERIES 

A. Graphical Gauge Model 

The determinant of a matrix was the key object discussed in IB]. Thus, we naturally start a technical description 
in this second paper with stating a new model in terms of determinants. 

Consider a square matrix, H , with elements Hah, a,b — I, ■ ■ ■ , N, and define a set of transformed (twisted) matrices, 
H{(t), determined by a set of fields, cr = {aat — ±l|a ^ b; a,b = 1, ■ ■ ■ ,N), hereafter referred to as gauge fields, 
according to the following rule: Hab{<y) = (JabHab and Haaicr) = Haa- The generalized ^-function of the matrix 
(understood as a generating function for Grassman variables correlations) is defined simply as the determinant, 

C(<T)^det(H(<T)). (1) 

Note that the Ihara function of a graph depends on a spectral parameter represented by a complex number lu. We 
add an additional set of binary spectral parameters, represented by the gauge field components and set to = 0. 

The matrix H also defines an undirected graph Q{H) with N nodes a & Qq. The nodes a and b are connected by an 
edge a = {a, b} e Qi when Hab 7^ or Hba 7^ 0. In other words the nodes a ^ Qq represent the diagonal elements Haa, 
whereas the edge {a, 6} G Qi corresponds to the pair of the off diagonal matrix elements Hab and Hab where at least one 
of the elements is nonzero. Hereafter we will also use a convenient notation a ~ 6 for {a, b} ^ Qi. Note that for defining 
the C-function we need only those components aab that are related to the edges of Q, i.e., a ^ b. Therefore, hereafter 
the gauge fields will include the relevant components, only, i.e. cr — (aab = ±l|a ^6; a ^ b; a,b — 1, ■ ■ ■ , N). 

The last comment allows the configurations cr to be interpreted as discrete gauge fields with the gauge group Z2 that 
reside on the graph Q in the full accordance with the terminology of the Lattice Gauge Theories [Hill [5], originally 
developed in the context of regular lattices. Recall that in Lattice Gauge Theories the components aab of a gauge field 
reside on the edges of the lattice and take values in the gauge group; the latter means that in our case the gauge group 
is Z2. It is important to note that generally aab '^ba, they rather satisfy the constraints ajja^^ab — 1 (with 1 being 
naturally the gauge group unit element). However, in our special case of the Z2 gauge group the constraints imply 
c^ha = Cab, and we can interpret the gauge field components as residing on the graph unoriented edges {a, &} = {6, a}, 
not on plaquettes as common in standard Lattice Gauge Theories considered on surface graphs ^. 

The Z2 gauge theory associated with matrix H and the graph Q{H), respectively, is stated simply as an average/sum 
of the gauge-field-dependent determinant over all possible configurations of the gauge field on the graph. The partition 
function of the model becomes 

2-1^^1 ^ C(t)^ / VcTdet{H{a)), (2) 

where the "integral" over the set of nonzero edges on the rhs is simply a convenient notation for the sum over 2l^il 
possible states of the discrete gauge fields and \Gi \ stands for the cardinality of Gi (i.e., number of edges of G)- 

Obviously one can think of any determinant on the rhs of Eq. ^ as of the one derived in the result of averag- 
ing/integration over Grassman variables 6, associated with the vertexes of G- Adopting the notation introduced in 
the first paper [B] (see also [12j ) we can recast the partition function ([2]) in a form 

Z = J Vcr J veveexp{SQ{e,e;(T)) , (3) 

So{9,9;<t) — ''^^ HaaOaOa+ ^ <y ab{HabO aOb + H^aObO a) ■ (4) 



^ A general graph does not have a notion of plaquettes, and therefore the gauge field curvature (intensity) that resides on the plaquettes 
may not be introduced. However, equipping the graph with an additional structure, namely a cyclic ordering of the edges attached 
to a vertex for all vertices a S 5oi which turns the graph into a so-called fat graph 1101 . allows to interpret Q as the 1-skeleton of a 
2-dimensional CW-complex that represents a Riemann surface, where the CW-complex 1111 is a space that can be obtained step-by- 
step via attaching cells of higher dimension. The set of points constitutes its 0-dimensional skeleton. Attaching 1-dimensional cells 
represented by the edges results in an un-oriented graph that constitutes the 1-dimensional skeleton. Attaching 2-dimensional cells 
represented by plaquettes results in the 2-dimensional skeleton. In our case the latter reproduces a Riemann surface and no more cells 
are attached. The gauge field intensity that resides in the 2-dimensional cells of the obtained CW-complex plays an important role in 
this case. 
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Following the terminology commonly accepted in the field theory and mathematical physics we call S'o the action of 
the Graphical Gauge Model. Note that, since the action in Eq. Q depends on the gauge field <t, it describes free 
fermions, interacting with the gauge field. The action of the pure gauge field in this model is zero. 



B. Monomer-Dimer Model 



The integrations/summations on the rhs of Eq. ^ obviously commutes, thus exchanging the order of integration, 
expanding vertex terms of the integrand in the series, utilizing the anti-commuting features of the Grassman variables 
and, finally, integrating over the binary gauge variables, one derives 

/" I?^e^«(^'«^-) = [] (l+wA^,) l[ {l+WabdaeaW, (5) 
•^^1 aeQo {a,b}egi 

where Wab = —HabHba and Wa = Haa- Expanding Eq. ([5| into a polynomial and integrating the resulting expression 
over the Grassman variables we find that only terms associated with valid monomer-dimer configurations survive (are 
nonzero), i.e. 

(6) 

where the set of tt consists of two sub-sets of binary 0, 1 variables defined on the vertexes of the graph and on the 
edges of the graph, respectively: tt = tt^, U TTe, tt^ = (tTq = 0, 1; a S Qq), and TTe = {iTab = 0, 1; {a, 5} € Qi). The last 
term on the rhs of Eq. ^ stated in terms of the Kroneker symbols describes the set of the monomer-dimer exclusions. 
In other words, a monomer-dimer configuration corresponds to a coloring of the graph (its vertexes and edges) in 
such a way that either at least one edge adjusted to the vertex is colored and then the vertex is not colored, or the 
adjusted vertexes are all uncolored and then the vertex is colored. 

Note that Eq. ([s]) after some obvious modification can be viewed as a definition of an effective action S{6, 6) that 
depends on the fermion variables only 




L 



D^e^o(M;-) ^gS(M)^ 5(^,0) = ^ U;,Ma+ ^abOaOaObeb, (7) 

^1 aeQo {a,b}eei 

As it usually happens in gauge theories, integration over the gauge field creates fermion interactions (second term in 
the action in Eq. ([t])). The interaction can be decoupled by introducing a Hubbard-Stratonovich field represented 
by another Z2 gauge field coupled to 9a0b and 9a9b- This results ina representation of the partition function of the 
Monomer-Dimer model in a form of an integral (sum) over the gauge field, with the integrand represented as a product 
of two gauge-field dependent Pfaffians. This representation will be studied in detail in the next paper of the series, 
with the focus on its applications to fat graphs. 



C. Oriented Disjoint cycle (Determinant) Series 

We further represent the integrand of the GGM partition function ^ in the following simple form 

^So(e,e;.) ^ Yl e--'-'- J] (e^-«""«^+^-«"^«" + (a,, - + HbaOb0a)) , (8) 

aeSo {Q,b}eei 

using straightforwardly the Grassman variables anticoagulation relations. Direct integration of Eq. ([s]) over the gauge 
variables implies 

p^e^o(M;<.) ^ -Q -Q (e^-^""^^+^-«"^''<. - (i/,,M. + i?.a^".0a)) ■ (9) 

aeGo {a,6}eei 

We further note that Eq. ^ can be represented as a sum of monoms in elements of H. Let us consider a monom 
which contains an off-diagonal element Hab but not its conjugate, Hab- Then, it is obvious (from the rules of the 
Grassman integration) that such a monom can only be associated with a directed disjoint cycle which contains the 
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directed segment {a,b), i.e. the monom should contain a product of the off-diagonal elements along the cycle and 
do not contain any of the respective conjugates. Moreover the product of the off-diagonal elements of H along the 
oriented disjoint cycle originates primarily from the expansion of the second product in Eq. ([9| in the series. Therefore, 
one concludes that Eq. (|9| can be represented as 



Z= J2 ^(^)' r{C)=aiC)dct{H\g\c) n (-^-b), 
CGODC(e) (a,b)ec 



aiC) 




n {^"^ 



{HabOadb + Hbadbda) , 



(a,6)GC 



(10) 



(11) 



H=0 



where |g\c denotes the restriction of iJ to ^ \ C. For a subgraph C C we denote by Q\C the maximal subgraph 
of G that has an empty intersection with C. Stated differently, Q \ C is represented by those edges of G that do not 
have common vertices with C. In Eq. ( 10 1 the det-term corresponds to direct integration over variables that do not 
belong the oriented disjoint cycle C. In essence, a{C) is a combinatorial factor which is calculated by straightforward 
counting. Expanding the integrand in Eq. (11) into a series over the square-bracket terms. One finds, that there are 

^ \ contributions associated with a product of k square-bracket terms along the oriented disjoint cycle, where 

|C| stands for the length of the oriented disjoint cycle measured in the number of segments/edges and 1 < fc < |C|, 
and each of them contributes (—1)'^"*"^ into a{C). Summing up all the nonzero contributions one derives. 



fe+i 



\c\ 
k 



= 1. 



(12) 



Substituting Eq. (12) into Eq. (10) we arrive at the desired expansion of the MD model partition function with the 

(13) 



coefficients represented by determinants 

Zmd = 



ceODC(e) 



r{C) ^ det{H\g\c) H (-^-"^ 

(a,b)eC 



An example of a family of oriented disjoint cycles for a sample graph is shown in Fig. [2] 



II. DETERMINANT AS A SERIES OVER MONOMER-DIMER CONTRIBUTIONS 

Eqs. (7-10,12) of 6^ represent the starting point for discussion of this Section. However, instead of following the 
path discussed in [6J, we make another non-BP choice of the gauge. 

The special gauges we will be using arc associated with the graph orientations d E 0{G), where 0{G) denotes 
the set of graph orientations. An orientation of G associates a direction ("arrow") on each edge, i.e., it is a pair of 
maps dj : Gi ^ Go with j = 0, 1 so that the edge a connects do{a) and di{a). For each edge there are two possible 
orientations: do{{a,b}) — a, di{{a,b}) = b and do{{a,b}) — b, di{{a,b}) = a. In particular ca,rd{0{G)) = 2^^^^^, 
where Nk — ca,rd{Gk) with A; = 0, 1 represent the number of edges and nodes. Therefore, orientation can be viewed 
as a binary variable that resides on the graph edges. The gauge associated with an orientation d £ 0{G), which is 
totally determined by specifying the set {7a6(<?)}a~fc of numbers that characterize the local ground states, is given by 

lab{d) (Haby^ for do{{a,b}), ^ab{d) = -{HabY^ for a = 9i({a, 6}), 

Hab = 1, Cab = 1/2, "f'ab = ^lab: Cab = {Hab)^^ (14) 

which simply means that we choose ^ab = '^{Hab)^^ depending on orientation, and the signs in front of Hab and iJ^a 
are always opposite. The rest of the parameters are determined by Eq. (12) of [B]. Note that the set of parameters 



{7o&(£^)}a~fc for a gauge choice given by Eq. (14) satisfy all the necessary requirements represented by Eq. (11) of 
[6]. Also note that two graph orientations d and d' are also related via a set cr of edge binary variables: we define 
'^{afi} = 1 if d{{a, 6}) = 9'({a, 6}), and cr{a,6} — ^li otherwise. In particular, a choice of some base graph orientation 
allows the graph orientations to be described using the edge binary variables cr. However, a generic unoriented graph 
is not equipped with a preferred choice of orientation. 



<i.p <i.p '^.p i> 




FIG. 2: Example of a graph (top) and its respective set of oriented disjoint cycles, consisting of 12 configurations. Oriented 
Disjoint cycles are shown in red. 



For any choice of a special gauge (14 1 the Grassman integral representation (Eq. (7) of ^) for the determinant of 
H can be represented in the following form 



det{H) = \ H {-HabHta) \ (l[Haa\ ^X^X J] ^ 

( \ 

\ah\ah\baX,ha \.ahX.ha \haX.ah 



X n 

{a.b}egi 



1 



HabHb, 



Hab 



ha 



(15) 



Veven states at {a, 6} odd states at (a, 6}/ 



which is explicitly independent on the choice of a special gauge. We further partition each factor on the second line 
of Eq. (151 that corresponds to an edge {a, 6} into a sum of two terms, referred to even and, odd according to the 
terminology introduced (and explained) in [B . Our next step is to expand the product of edge terms in the integrand 
of Eq. ( 15 ) into a polynomial over the odd states, followed by performing integration over the edge Grassman variables 



that correspond to the odd contributions. For a given choice of the local odd states we denote hy C d Q the subgraph 
of Q formed by the edges, where the odd terms [the third or forth term in the second line of Eq. ( 15 )] have been chosen. 



We start with demonstrating that all vertices of the subgraph C C Q have the valence two, i.e., C is represented by a 
disjoint union of simple loops. This follows from the fact that the expression in the exponent in Eq. (15 1 is actually 
a product of two linear combinations of the original Grassman variables and, therefore 



(Ha 



:i.Tio ^^'^ = 1 + (Haa)-' ^ba ^b'a. 

beGo b'eGo 



(16) 



Consider a vertex a G Cq- For the integral over the local vertex variables dXadXa not to vanish the integrand should 
contain each of the local variables Xba and Xba with b ^ a exactly once. These local variables in the integrand originate 
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from the odd terms, described above, from the even terms [second contribution in the second Hne of Eq. ( 15 )] and from 
the relevant exponential terms represented by Eq. (16 1. If an edge {b,a} € Gi also belongs to C we have either the 
odd term {Hab)~^XabXba or {IIi,a)~^XbaXab in the integrand. Consider the first option (the second option is considered 
in a similar way), the local conjugate variable Xba can originate only from the vertex term given by Eq. (16 1 and 



is represented by a contribution (Haa) XbaXh'a- The variable Xfa conjugate to the the variable Xb'a can originate 
only from an odd edge contribution, namely {Hb'a)^^Xb'aXb'a, which implies that {b',a} G Ci. The edges {b,a} and 



{b',a} are the only edges adjacent to the node a that belong to Ci since the vertex term [Eq. (16)] that provides 
the conjugate, local variables contains products of only two Grassmans. Consideration of the other odd edge term 
{Hi,a)~^XbaXab leads to a similar result, but with the opposite orientation. Therefore, any vertex a £ C has the 
valence two, and, therefore C S OT)C{Q). 



Consider a simple oriented loop (oi , 



i), where naturally aj ^ Oj+i and a„ ^ ai. The associated contribution 



given by the integral of the loop edge variables of the product of the edge and vertex contributions has a form 



/n 



n 

X J_ J_ Xaj-f-iaj Xaj^iaj 



n „ 

Xoj_iQ,jXoj^l(2jXQj + lCljXoj_iC[.j 



(17) 



and in Eq. (17) we use a cyclic convention j + n ~ j. The first equality in Eq. (17) is obtained by performing 



permutations in the following way. We start with moving the Grassman Xa4,a3 in the integrand by two places to the 
left, followed by moving the combination Xa2a3Xaia3 to combine it with the combination Xa4a3Xo2a3 in the product 
over J, which corresponds to the value j = 3; after that we permute the Grassmans Xa3a2 and Xaaa^- The overall 
permutation provides a ( — 1) sign factor. Repeating a similar operation (n — 1) times (including the first explicitly 
described operation) results in the first equality. The second equality follows from the fact that each of n Grassman 
integrals in the intermediate expression is equal to (—1). 

Due to Eq. (17) the resulting expression for the determinant adopts a form 



det{H) 



E (- 

ceODC(e) 



T) 



dog(C) 




(18) 



Zl{g\C)=[ Yl i-HabHta)\ i n 

Ja,6}e(ec)l / KaeiGch 



h' r^a 



T^xGcPxQc n i+(^oa)"' ^^'A n 



XahXabXbaXba 



(19) 



where deg(C) denotes the number of connected components in C . Stated differently, an element C G ODC(C/) is 
represented by a disjoint union of oriented simple loops with deg(C) denoting the number of simple loops in C . In 
Eq. (18) XQc^ XGc denote the edge Grassman variables restricted to the subgraph Qc, formed by the edges of Q that 
do not belong to C. In deriving Eq. ( [T9| we have also made use of Eq. (16) to replace the exponential vertex terms 
with their polynomial counterparts. 



It is now straightforward to check [by expanding the integrand in Eq. ( 19 1 into a polynomial followed by performing 
integration over the Grassman variables in Eq. (19)] that Zi{Q') is nothing else than the partition function (|6| of the 
monomer-dimer model on the graph Q. Consider an edge {a, 6} g {Gc)i- The Grassman variables Xab, Xba, Xab Xba 
whose product provides a nonzero contribution to the integral over the edge variables can originate from the vertex or 
edge terms in Eq. (19). If they originate from the edge term, then combining with the corresponding edge prefactor 



[from the first line in Eq. ( 19 1] we obtain the contribution 



HabHba / dXabdXabdXbadXba 



XabXabXbaXba 
HabHba 



= 1, 



(20) 
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if they come from the vertex terms associated with the vertices a and h, then combined with the corresponding vertex 
pref actors, the contribution has a form 



HaaHbb J dXabdXabdXbadXba{Hbb) XabXab{Haa) (21) 

We call such an edge a dimer. Obviously, any node can have not more than one dimer edge attached to it. The 
nodes that do not have dimers attached to them are referred to as monomers. A monomer node a does not provide 
the Grassmans associated with the vertex term and, therefore, the prefactor term Haa is not compensated. A dimer 
{a, 6} does not provide the edge terms and, therefore, the edge prefactor {—HabHba) is not compensated. It is easy to 
see that any configuration of monomers and dimers that provides a non-zero contribution to the Grassman integral 



in Eq. (19 1 satisfies the monomer-dimer matching rules. Therefore, Zi{Q \ C) represents the partition function of the 
monomcr-dimer model with the monomer and dimer weights Wa — Haa and Wab — —HabHba, respectively. 
Summarizing, 

vc, g' = g\C: Zi{g') = ZuBiG')- (22) 

which implies 



det(iJ)- ^ r{C), r(C) = (-l)'^^s(C) "Q (^Hab)ZMB{g\C) (23) 
ceODC{g) (a,fc)ec 

To conclude, we just showed that the determinant of a matrix can be represented in terms of a series over disjoint 
oriented cycles of the underlying graph, with each term of the expansion being proportional to the partition function 
of the monomer-dimer model defined on the remainder of the graph, i.e., after the cycles, as well as all edges connected 
to their vertices are removed. 



Comparing Eq. (13 1 with Eq. (23 1 one finds that in a sense one is an inverse of the other. While the former 
expresses the partition function of the MD model on the graph in terms of an expansion over the determinants (each 
corresponds to a directed disjoint cycle), the later does exactly the opposite by expressing the determinant as a series 
over the partition functions of the MD models each associated with the exclusion of a directed disjoint cycle. More 
details on this relation are given in Appendix \K\ 

We complete this Section by addressing the issue of the gauge invariance of the simple-loop decomposition. To that 



end we twist the matrix H as described at the beginning of Section I A i.e. introducing the matrix H{cr), twisted 
by the gauge field cr as Habic) — (JabHab ior a ^ b and Haa{(^) — Haa- Applying Eq. (23) to H{cr), recaUing the 
definition of the ^-function ([T]), and noting that the partition functions Zmd{G \ C) are obviously invariant with 
respect to the twisting we obtain the following decomposition for the ^-function: 

C(^) = det(i7(*T)) = (-1)"=""^^ n {-<^abHab)ZMT,{g\C)^ ^ r{C) J] ^ab (24) 

C&ODC(G) {a,b)eC CeODC(g) (a,b)eC 

Therefore, r(C) can be viewed as the coefficients in the expansion of the ^-function C(''') in the gauge field cr and, 
therefore they do not depend on a particular way they are evaluated. 



III. SUMMARY AND CONCLUSIONS 



To summarize, this manuscript reports new relations between the partition function Zmd{Q) of the monomer- 
dimer model, defined on an arbitrary graph Q and the corresponding determinant of the matrix H and its minors, 
constructed from the monomer-dimer weights on the graph. We have formulated a Graphical Gauge Model (GGM) on 
a graph, stated in terms of Grassman variables and binary gauge fields, so that all the relations reported in the paper 
follow in a straightforward way via simple and natural manipulations (reparametrizations and integrations) over the 
partition function of the GGM. Some results of this paper are also linked to the discussions in the first paper of the 
series [5]. In particular, we show here that the expression for a determinant as an expansion over directed disjoint 
cycles is related to the Loop Series approach of The difference comes from different gauge choices. 

In spite of the progress in understanding relations between determinants, loops and matchings (i.e. valid configura- 
tions of the monomer-dimer problems), there are still many important challenges left for future analysis. We conclude 
with mentioning some of these "natural" challenges. 
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• Given the prominent role the determinants play in the classical studies of the dimer models on planar graphs and 
graphs embedded in Riemann surfaces of finite genus [13l [TH [151 [El Ull UHl ES] j one suggests that it should be 
important to analyze the consequences of the monomer-dimer, determinant, loops and GGM relations discussed 
above for planar and surface graphs, also extending the results of [ST. 

• All the Loop Series related constructions for graphical models, introduced so far in [H [51 UHl [S] and this 
manuscript, express the partition functions as series over sub-graphs. On the other hand, the well-known formula 
Indet(if) = Trln(77), and related famous expressions for the log-partition function of the Ising model on a 
planar graph [22 , suggests that a multiplicative expansion that represents the partition functions as a product 
over sub-graphs, may also exist, at least for some class of graphical models. Exploring possible multiplicative 
decompositions constitutes an important theoretical and algorithmic challenge. 

• One general technical conclusion of the paper is related to the use of Berezin integrals |12j . Our approach shows 
that the Grassman-integration technique can be useful for deriving quantitative exact relations in graphical 
statistical problems of computer science, operation research, and information theory. Obviously, the two papers 
of the series present only the first step in this direction. A possible extension of this approach, worth a future 
exploration, would be to develop a more general super-symmetrical and cr-models based approach, in the spirit 
of |23| . combining normal and Grassman integrations. 

• To a large extent, the practical utility of the determinant and cycle series discussed in the paper is yet to be 
determined. In particular, it remains to be seen weather the reported cycle series allows an efhcient deterministic 
approximation for the monomer-dimer model partition function. We speculate that an algorithmic extension of 
our results may lead to the development of novel Fully Polynomial-Time Approximation Schemes (FPTAS) for 
various hard, #P, weighted counting problems (see e.g. j24j for a sample FPTAS example discussed recently). 
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APPENDIX A: CORRESPONDENCE BETWEEN MONOMER-DIMER MODEL AND DETERMINANT 

SERIES 



In this Appendix we establish relation between Eq. ( 13 1 and Eq. ( 23 1 in a somehow straightforward way. 



A particular strength of the decomposition of the determinant Eq. (23 1 is its naturality, i.e., it is valid for any graph 



g associated with some matrix H . In particular it can be written for any subgraph Q' <Z G . To see the advantages 
of naturality in a more clear way we introduce the following notation Xc — det{Q \ C) and Yc — ^md(^/ \ C) where 
C £ ODC{Q) is an (oriented) simple loop in Q. For our purposes it is also convenient to introduce an oriented graph 
C.{Q) of simple loops, whose nodes are simple loops C G ODC{Q), i.e., Cq{Q) = ODC{Q). The set of links of an 
oriented graph Ci{Q) C Cq{Q) x Co{Q) is naturally a subset Ci{Q) C ODC(Q) and is defined as follows. We say that 
(C, C") is an oriented link (a connecting arrow goes from C" to C), i.e. (C, C") G Ci{Q), if C C C". 

The reason why the graph C. =C.{Q) has been introduced is that C. is the oriented graph associated with the linear 
relation (matrix) that expresses the set {Xc}c£Co of partition functions in terms of the set {IcIceCo of partition 
functions. To see that we recast Eq. ( [23| for an arbitrary subgraph Q' d Q 

det{g')= J2 (-1)''''*''' n i-Hab)ZMBig'\C). (Al) 
ceonc(e') (a,fc)ec 



Applying Eq. ( Al I for all Q' = Q \ C and making use of the introduced notation we arrive at 



(c,c')eCi 

C {a,b)eC'\C 



(A2) 



C'DC 



Note that Rcc" 7^ 0, if and only if (C, C") G Ci, i.e., (C, C") is an edge of the oriented graph C, which means that C. 
is the oriented graph associated with the ma trix Rcc- 

The oriented disjoint cycle expansion ( 13 1 for Zmd is obtained by expressing the inverse matrix R^^, as a sum 
over the oriented (i.e., orientation on the path should be compatible with the orientation on the graph) paths on the 
associated graph C. : 



pePC. j=o 

Po=C',p,(p)=C 

= n (-Hab) E (-l)'(rt+dcg(C')-dog(C) ^ -Q (_^^^), (A3) 

{a,b)eC'\C pePC. {aM}eC'\C 



In deriving Eq. ( A3 ) we made use of the fact Rcc ~ 1 for the diagonal elements and the expression for the off-diagonal 
components ( A2 1 . In particular the specific form of Rcc implies that the contributions of different paths are the 



same up to a sign. The last equality in Eq. (A3 1 is obtained by an explicit computation of the combinatorial factor 



dcg(C')-dog(C) 

E (- 



.^)/(p)+dcg(c')-dcg(c)_^ deg(C') - deg(C)) 



(A4) 



1=0 



where Af{l', N) is the number of ways one can put N objects into / boxes with each box containing at least one object. 
Applying Eq. (A4| to C = and recalling the meaning of the notation Xc and Yq we arrive at Eq. (13). 
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APPENDIX B: EXPANSION OF A DETERMINANT AND SUMMATION OVER THE GAUGES 



In this Appendix we present an alternative derivation of the decomposition ( 23 1 of a determinant into a sum over 



the oriented disjoint cycles with the individual contributions expressed in terms of the partition functions of the 
Monomer-Dimer (MD) models defined on the proper subgraphs of G- 

First of all we note that the loop decomposition (Eq. (22) of [5]) is valid in any gauge, i.e., for any choice of the 
set {7a6}a~fc, provided that the summation over generalized loops C e GL{Q) is extended to the summation over all 
subgraphs Q' G Q. In the case of a BP gauge the latter summation is restricted to the summation over the generalized 
loops, since the BP gauge ensures the vanishing of the rest of the contributions. Multiplying the relative contributions 
r(C, C") with the prefactor and changing the order of the summations we recast the loop decomposition in a form 



ceODC(e) 



r{C) 



cce' 
e'ce 



(Bl) 



of a decomposition in oriented disjoint cycles. Note that in this notation the BP contribution corresponds to the 
empty simple loop and empty subgraph. Note that strictly speaking the loop series depends on the gauge choice. 
However, the gauge freedom (among the gauges we are dealing with) belongs to the boson subspace, which implies 
that the coefficients r{C) in Eq. ( |B1|) sh ould be gauge invariant. This issue is addressed at the end of section [ll| 

We will consider the loop series (Bl) for all 2^^ special gauges associated with the graph orientations (they are 
given by Eq. (14l) and average it with an equal weight of 2"^^. This is a legitimate procedure since the sum of all 
terms in a loop series is naturally gauge invariant. We also note that for given C S ODC{Q) a particular choice of a 
subgraph C G Q' G G can be described by a particular configuration of a set cr S Mc of binary variables that reside 
on those edges of Q that do not belong to C. Namely, ctq = —1 for a & Q' (painted edge that correspond to local 
even excited state) and ctq, = 1 otherwise (local ground state). Combining these arguments with the expressions for 
the ingredients of the loop expansion )Eqs. (23) and (24) in |6]) we arrive at 



cce' 

Z = 2-^^(e) E E E Z{Q',C;d) 

9eo{g) ceODC{g) S'ce 

= ^ EE (-l)'^°^^'^^2"(^i(^)-^i(<^)) 

C£ODC(g)deo{g)TeMc (c, 
= 2-Ni{g) ^ ^__-[^-jdog(c)2-(^Vi(e)-7Vi(c)) 

ceODC{g) 



n ^-'^ n \-Haa-Y,^lba{d)r^^ba\ 
,d)eC ae(e\C)o V b~a ) 



n ^--i E E n [-iiaa-Y.^lUd{{aM))r^^ba\ 
(c4)ec o-eMc aeo(g) ae(g\c)o \ b^a ) 



(B2) 



Comparing Eq. (B2| with Eqs. (Bl) we see that the decomposition in simple loops (23) is reproduced if we define 



^MD(e\C) = 2-(^i(e)-^i(^)) E 2""^'^^^ 



E n 

seo(e) ae(e\c)o 



— Una — 



(76a(a({a,6})))- 



(B3) 



The only thing we need to show at this point is that the expression in Eq. ( B3 1 reproduces the partition function of 



the MD model on the graph G\C. This is achieved by performing the summation over the binary variables d G 0{G)- 
The desired result follows from an obvious property 



E {lba{d{{a,b}W 

9{{aM) 



Cba = 0, 



E {lba{d{{a,bm-' {j,,{d{{aM))y 

9{{aM) 



f^baf^ab 



-HnhH] 



ab^ ba ; 



(B4) 



where both sums in Eq. (B4| contain two terms that correspond to two possible values of the orientation c?({a, 6}) of 
the edge {a, b}. A choice of a diagonal term in the parenthesis in Eq. ( |B3[ ) corresponds to having a monomer on the 
node a with the weight Haa- It follows from Eq. (B4| that the off-diagonal terms should always go in pairs, each pair 
(7a67ba)~^ corrcspouds to having a dimer on the link {a, 6}, whose weight is 



-HabHba- It also follows from Eq. (|B4| 
that the sum over orientations in Eq. (B3l does not depend on cr and, therefore, the sum over cr just cancels out the 
first prefactor in Eq. ( B3 1 . This completes the proof. 



